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EDITORS’ PREFACE. 


Tue volume called Higher Mathematics, the first edition 
of which was published in 1896, contained eleven chapters by 
eleven authors, each chapter being independent of the others, 
but all supposing the reader to have at least a mathematical 
training equivalent to that given in classical and engineering 
colleges. ‘The publication of that volume is now discontinued 
and the chapters are issued in separate form. In these reissues 
it will generally be found that the monographs are enlarged 
by additional articles or appendices which either amplify the 
former presentation or record recent advances. This plan of 
publication has been arranged in order to meet the demand of 
teachers and the convenience of classes, but it is also thought 
that it may prove advantageous to readers in special lines of 
mathematical literature. 

It is the intention of the publishers and editors to add other 
monographs to the series from time to time, if the call for the 
same seems to warrant it. Among the topics which are under 
consideration are those of elliptic functions, the theory of num- 
bers, the group theory, the calculus of variations, and non- 
Euclidean geometry; possibly also monographs on branches of 
astronomy, mechanics, and mathematical physics may be included. 
It is the hope of the editors that this form of publication may 
tend to promote mathematical study and research over a wider 
field than that which the former volume has occupied. 


December, 1905. 
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AUTHORS) PREEAGE. 


In republishing this short treatise in book form the author 
solicits criticism but offers no apology. The type of the book 
he has sought to imitate is that shown in the ‘“ mathematical 
tracts” of the late Sir George B. Airy. The brevity and the con- 
crete illustrations of these “ tracts ”” have Served very effectively 
in introducing students to a number of the more difficult fields 
of applied mathematics; and it is hoped that this treatise will 
serve a similar end. 

The theory of probability and the theory of errors now con- 
stitute a formidable body of knowledge of great mathematical 
interest and of great practical importance. Though developed 
largely through applications to the more precise sciences of as- 
tronomy, geodesy, and physics, their range of applicability extends 
to all of the sciences; and they are plainly destined to play an 
increasingly important réle in the development and in the appli- 
cations of the sciences of the future. Hence their study is not 
only a commendable element in a liberal education, but some 
knowledge of them is essential to a correct understanding of 
daily events. 

No special novelty of presentation is claimed for this work; 
but the reader may find it advantageous to know that a definite 
plan has been followed. This plan consists in presenting each 
principle, first, by means of a simple, concrete example; passing, 
secondly, to a general statement by means of a formula; and, 
thirdly, illustrating applications of the: formula by concrete 
examples. Great pains have been taken also to secure clear and 
correct statements of fundamental facts. If these latter are 
duly understood, the student needs little additional aid; if 
they are not duly understood, no amount of aid will forward him. 

The passage from the elementary concrete to the advanced 
abstract may appear to be abrupt to the reader in some cases. 
It is hoped, however, that any large gaps may be easily bridged 
and that any serious difficulties may be easily overcome by means 
of the references given to the literature of the subject. In any 
event the student will find that in this, as in all of the more ardu- 
ous sciences, his greatest pleasure and his highest discipline will 
come from bridging such gaps and from surmounting such diffi- 
culties. 

WasHINGTON, D. C., December, 1905. 
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PROBABILITY AND THEORY OF ERRORS. 


ART. 1. INTRODUCTION. 


IT is a curious circumstance that a science so profoundly 
mathematical as the theory of probability should have origi- 
nated in the games of chance which occupy the thoughtless and 
the profligate.* That such is the case is sufficiently attested 
by the fact that much of the terminology of the science and 
many of its familiar illustrations are drawn directly from the 
vocabulary and the paraphernalia of the gambler and the trick- 
ster. It is somewhat surprising, also, considering the antiquity 
of games of chance, that formal reasoning on the simpler 
questions in probability did not begin before the time of Pascal 
and Fermat. Pascal was led to consider the subject during the 
year 1654 through a problem proposed to him by the Chevalier 
de Méré, a reputed gambler.t The problem in question is 
known as the problem of points and may be stated as follows: 
two players need each a given number of points to win ata 
certain stage of their game; if they stop at this stage, how should 
the stakes be divided? Pascal corresponded with his friend 
Fermat on this question ; and it appears that the letters which 
passed between them contained the earliest distinct formulation 
of principles falling within the theory of probability. These 


* The historical facts referred to in this article are drawn mostly from Tod- 
hunter’s History of the Mathematical Theory of Probability from the time of 
Pascal to that of Laplace (Cambridge and London, 1865). 

+‘ Un probléme relatif aux jeux de hasard, proposé aun austére janséniste 
par un homme du monde, a été l’origine du calcul des probabilités.” Poisson, 
Recherches sur la Probabilité des Jugements (Paris, 1837). 
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acute thinkers, however, accomplished little more than a correct 
start in the science. Each seemed to rest content at the time 
with the approbation of the other. Pascal soon renounced 
such mundane studies altogether; Fermat had only the scant 
leisure of a life busy with affairs to devote to mathematics; 
and both died soon after the epoch in question,—Pascal in 
1662, and Fermat in 1665. 

A subject which had attracted the attention of such dis- 
tinguished mathematicians could not fail to excite the interest 
of their contemporaries and successors. Amongst the former 
Huygens is the most noted. He has the honor of publishing 
the first treatise* on the subject. It contains only fourteen 
propositions and is devoted entirely to games of chance; but it 
gave the best account of the theory down to the beginning of 
the eighteenth century, when it was superseded by the more elab- 
orate works of James Bernoulli,t Montmort,t and De Moivre.§ 
Through the labors of the latter authors the mathematical 
theory of probability was greatly extended. They attacked, 
quite successfully in the main, the most difficult problems; 
and great credit is due them for the energy and ability dis- 
played in developing a science which seemed at the time to 
have no higher aim than intellectual diversion. Their names, 
undoubtedly, with one exception, that of Laplace, are the most 
important in the history of probability. 

Since the beginning of the eighteenth century almost every 
mathematician of note has been a contributor to or an expos- 
itor of the theory of probability. Nicolas, Daniel, and John 
Bernoulli, Simpson, Euler, d’Alembert, Bayes, Lagrange, Lam- 
bert, Condorcet, and Laplace are the principal names which 
figure in the history of the subject during the hundred years 

* De Ratiociniis in Ludo Alew, 1657. 

+ Ars Conjectandi, 1713. 

} Essai d’Analyse sur les Jeux de Hazards, 1708. 

§ The Doctrine of Chances, 1718. 

| Todhunter says of Montmort, for example, “In 1708 he published his 


work on Chances, where with the courage of Columbus he revealed anew world 
to Mathematicians.” 
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ending with the first quarter of the nineteenth century. Of 
contributions from this brilliant array of mathematical talent, 
the Théorie Analytique des Probabilités of Laplace is by far 
the most profound and comprehensive. It is, like his Mé- 
canique Céleste in dynamical astronomy, still the most elabo- 
rate treatise on the subject. An idea of the grand scale of the 
work in its present form* may be gained by the facts that the 
non-mathematical introductiont covers about one hundred and 
fifty quarto pages; and that,.in spite of the extraordinary 
brevity of mathematical language, the pure theory dnd its ac- 
cessories and applications require about six hundred and fifty 
pages. 

From the epoch of Laplace down to the present time the 
extensions of the science have been most noteworthy in the 
fields of practical applications, as in the adjustment of obser- 
vations, and in problems of insurance, statistics, etc. Amongst 
the most important of the pioneers in these fields should 
be mentioned Poisson, Gauss, Bessel, and De Morgan. Nu- 
merous authors, also, have done much to simplify one or an- 
other branch of the subject and thus bring it within the range 
of elementary presentation. The fundamental principles of 
the theory are, indeed, now accessible in the best text-books 
onalgebra: and there are many excellent treatises on the pure _ 
theory and its various applications. 

Of all the applications of the doctrine of probability none 
is of greater utility than the theory of errors. In astronomy, 
geodesy, physics, and chemistry, as in every science which at- 
tains precision in measuring, weighing, and computing, a 
knowledge of the theory of errors is indispensable. By the aid 
of this theory the exact sciences have made great progress dur- 


* The form of the third edition published in 1820, and of Vol. VII of the 
complete works of Laplace recently republished under the auspices of the 
Académie des Sciences by Gauthier-Villars. This Vol. VII bears the date 1886. 

+ ** Cette Introduction,” writes Laplace, ‘‘est le developpement d'une Legon 
sur les Probabilités, que je donnai en 1795, aux Ecoles Normales, ow je fus ap- 
pelé comme professeur de Mathématiques avec Lagrange, par un décret de la 


Convention nationale.” 
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ing the nineteenth century, not only in the actual determination 
of the constants of nature, but also in the fixation of clear 
ideas as to the possibilities of future conquests in the same di- 
rection. Nothing, for example, is more satisfactory and in- 
structive in the history of science than the success with which 
the unique method of least squares has been applied to the 
problems presented by the earth and the other members of the 
solar system. So great, in fact, are the practical value and 
theoretical importance of the method of least squares, that it is 
frequently mistaken for the whole theory of errors, and is 
sometimes regarded as embodying the major part of the doc- 
trine of probability itself. 

As may be inferred from this brief sketch, the theory of 
probability and its more important applications now constitute 
an extensive body of mathematical principles and precepts. 
Obviously, therefore, it will be impossible within the limits of 
a single condensed monograph to do more than give an out- 
line of the salient features of the subject. It is hoped, how- 
ever, in accordance with the general plan of the volume, that 
such outline will prove suggestive and helpful to those who 
may come to the science for the first time, and also to those 
who, while somewhat familiar with the difficulties to be over- 
tome, have not acquired a working knowledge of the subject. 
Effort has been made especially to clear up the difficulties of 
the theory of errors by presenting a somewhat broader view of 
the elements of the subject than is found in the standard 
treatises, which confine attention almost exclusively to the 
method of’ least squares. This chapter stops short of that 
method, and seeks to supply those phases of the theory which 
are either notably lacking or notably erroneous in works 
hitherto published. It is believed, also, that the elements here 
presented are essential to an adequate understanding of the 
well-worked domain of least squares.* 


*The author has given a brief but comprehensive statement of the method 


of least squares in the volume of Geographical Tables published by the Smith- 
sonion Institution, 1894. 
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ART. 2. PERMUTATIONS. 


The formulas and results of the theory of permutations 
and combinations are often needed for the statement and so- 
lution of problems in probabilities. This theory is now to be 
found in most works on algebra, and it will therefore suffice 
here to state the principal formulas and illustrate their mean- 


ing by a few numerical examples. 


- The number of permutations of 2 things taken 7 in a group 


is expressed by the formula 


(2), = n(n — 1)n—2)...(a— 


7-1). (1) 


Thus, to illustrate, the number of ways the four letters a, 6, 
¢, @ can be arranged in groups of two is 4. 3 = 12, and the groups 


are 


ah, ba, ac, ca, ad., da, bc, cb, bd, db cd, de. 


Similarly, the formula gives for 
fe ean == 2, pk Ca Ce 
m= 7“ r=3; (7)s 7-6.5 


II 


cad vt 
= 210, 


n=10 “ r=6, (10),= 10.9.8.7.6.5 = 151200. 


The results which follow from equation (1) when # andr 


do not exceed 10 each are embodied in the 


VALUES OF PERMUTATIONS. 


following table: 


10 9 8 4) 6 | 5 Pe Ne Seg eB 
t 10 9 8 7 6 5| 4] 3)2) 12 
2 go 72 £6 42 30 20°} 732 6 he 
3 720 504 336 210 120 60 | 24 6 
4 5040 3024 1680 840 360:} 120°] 2 
5 30240 15120 6720°| 2520 720 | 120 
6 151200 60480 20160 5040 720 
yf 604800 | 181440 40320 5040 
8 1814400 | 362880 | 40320 
9 3628800 | 362880 
10 3628800 
a Sa — Maas 
| Sp | 9864100 986409 | 109600 | 13699 | 1956 | 325 | 64 | 15 | 4] I 


‘The use of this table is obvious. Thus, 


the number of per- 


mutations of eight things in groups of five each is found in the 
fifth line of the column headed with the number 8. It will be 
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noticed that the last two numbers in each column (excepting 
that headed with 1) are the same. This accords with the for- 
mula, which gives for the number of permutations of x things. 
in groups of x the same valué as for ~ things in groups of 
(1 —1). It will also be remarked that the last number in each 
column of the table ‘is the factorial, z!, of the number z at the 
head of the column. For example, in the column under 7, the 
last numbers: 5040 = "142. 3.45 Ose a 

The total number of permutations of # things taken singly, 
in groups of two, three, etc., is found by summing the numbers 
given by equation (1) for all values of x from 1 to” Calling 
this total or sum S,, it will be given by 


Sy 2%) (2) 
To illustrate, suppose z = 3, and, to fix the ideas, let the 
three things be the.three digits 1, 2, 3. Then from the above 


table it is seen that S, = 3 4+ 6-+4+6=15; or, that the number 


of numbers (all different) which can be formed from those dig- 
its is fifteen, These numbers are 1, 2, 3; 12, 13, 21, 23, 31, 92% 
Pod, UA FUR, Oak. Ghee Cae, 
The values of S$, for z= I, 2,,.... 10 are piven undereiie 
corresponding columns of the above table. But when z is 
‘large the direct summation indicated by (2) is tedious, if not 
impracticable. Hence a more convenient formula is desirable. 
To get this, observe that (1) may be written 
n\ 

(x), = wry)’ ; (1) 
if 7 is restricted to integer values between 1 and (z — 1), both 
inclusive. This suffices to give all terms which appear in the 
right-hand member of (2), since the number of permutations 
for y= (z— 1) is the same as forr =~. Hence it appears 
that 


! 
Fic 7 one ear 


ee ee 


a 


ee ee ee ee ee 8 
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But as x increases, the series by which #! is here multiplied 
approximates rapidly towards the base of natural logarithms; 
that is, towards 


é = 2.7182818 +, log e = 0.4342945. 
Hence for large values of x 
S, = n\e, approximately.* (3) 
To get an idea of the degree of approximation of (3), sup- 


pose z=g9. Then the computation runs thus (see values in 
the above table): 


log 
9! = 362880 5.5597630 
e 0.4342945 
g!e= 986410 = 5-9940575 
S, = 986409 by equation (2). 
The error in this case is thus seen to be only one unit, or 
about one-millionth of S;.+ 


Prob. r. Tabulate a list of the numbers of three figures each 
which can be formed from the first five digits 1,...5. Howmany 
numbers can be formed from the nine digits? 

Prob. 2. Is S$; always an odd number for 2 odd? Observe 
values of Sy in the table above. 


ART. 3. COMBINATIONS. 


In permutations attention is given to the order of arrange- 
ment of the things considered. In combinations no regard is 
paid to the order of arrangement. Thus, the permutations of 
the letters a, 6, c, din groups of three are 


(abc) (abd) bac bad ach (acd) cab cad 
adb ade dab dac bea (bcd) cha chad 
bda bade dba dbe cda cdb dea dadeb 


* See Art. 6 for a formula for computing #! when x is a large number. 

+ When large numbers are to be dealt with, equations (1)’ and (3) are easily 
managed by logarithms, especially if a table of values of log (x!) is avaiiable. 
Such tables are given to six places in De Morgan’s treatise on Probability in 
the Encyclopedia Metropolitana, and to five places in Shortrede’s Tables 


(Vol. I, 1849). 
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But if the order of arrangement is ignored all of these are 
seen to be repetitions of the groups enclosed in parentheses, 
namely, (abc), (abd), (acd), (bcd). Hence in this case out of 
twenty-four permutations there are only four combinations. 


A general formula for computing the number of combina- 
tions of z things taken in groups of 7 things is easily derived. 
For the number of permutations of z things in groups of 7 is 
by (1) of Art. 2 

(7), = n(x — 1)(#— 2)...(a—7r+1); 
and since each group of 7 things gives 1.2.3...7= 7! per- 
mutations, the number of combinations must be the quotient 
of (”), by r!. Denote this number by C(z),. Then the gen- 
eral formula is 


Cn = n(n — 1)(# —2)...(#—7-+ 1) 


7! 


(1) 


This formula gives, for example, in the case of the four let- 

ters a, 6, c,d taken in groups of three, as considered above, 

dees es 

C(4), = a ae 

Multiply both numerator and denominator of the right-hand 
member of (1) by (z — 7)! The result is 


4. 


a} 
al ater (x—7)! (1)! 
which shows that the number of combinations of z things in 
groups of 7 is the same as the number of combinations of 
things in groups of (z— 7). Thus, the number of combina- 
tions of the first ten letters a, 6,c...7 in groups of three or 


seven is 
10! 
317! = 120. 
The following table gives the values C(x), for all values of 
wz and 7 from fr to Io. 


The mode of using this table is evident. For example, the 
number of combinations of eight things in sets of five each is 
found on the fifth line of the column headed 8 to be 56. 


COMBINATIONS, 15 


VALUES OF COMBINATIONS. 


I 10 9 8 7 6 5 A r3|"2-1_5 
2 45 36 28 2Ey\- 55 |} tO Onegai ay 
3 120 84 56 Bou 20 a 4|1 

4 210 | 126 | -70 | 35 | 15 I 

5 252 | 126 56 21 6 I 

6 210 | 84| 28 al ies 

yi 120 | 36 8 I 

8 45 9 I 

9 10 I 

10 I 


Se 1023 | 511) 255 | 127 | 63.]| 32 | 15 | 74 3 


It will be observed that the numbers in any column show 
a maximum value when z is even and two equal maximum 
values when zis odd. That this should be so is easily seen 
from (1)’, which shows that C(z), will be a maximum for any 
value of x when r!(z — 7)! isa minimum. For x even this is 
a minimum for 7 = $7; while for z odd it has equal minimum 

- values for 7 = 4(z— 1) andr=}(z-+1). Thus, 

maximum of C(z),= fy for even, 
; (2) 
n! 

——— for # odd; 


m+ ards 


7 
hic ey 
2 


2 


= 


The total number of combinations of z things taken singly, 
in groups of two, three, etc., isfound by summing the numbers 
given by (1) for all values of 7 from 1 to # both inclusive. 
Calling this total or sum S., 

5, == 2O(R),. 
The same sum will also come from (1) by giving to 7 all values 
from 1 to (z — 1), both inclusive, summing the results, and in- 
creasing their aggregate by unity. Thus by either process 


+...-a2+1. 


n(n — 1), n(n — 1" 3) 


ee 12 1.2.3 
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The second member of this equation is evidently equal to 
(1+1)*—1. Hence 
S,= SC) eer is (3) 
The values of S, for values of 2 and x from I to 10 are given 
under the corresponding columns of the above table. 


Prob. 3. How many different squads of ten men each can be 
formed from a company of roo men ? 

Prob. 4 How many triangles are formed by six straight lines 
each of which intersects the other five ? 

Prob. 5. Examine this statement: “In dealing a pack of cards 
the number of hands, of thirteen cards each, which can be produced 
is 635 013559 600. But in whist four hands are simultaneously held, 
and the number of distinct deals . . . would require twenty-eight 
figures to express it.” * 

Prob. 6. Assuming combination always possible, and disregarding 
the question of proportions, find how many different substances 
could be produced by combining the seventy-three chemical ele- 
ments. 


ArT. 4. DIRECT PROBABILITIES. 


If it is known that one of two events must occur in any 
trial or instance, and that the first can occur in a ways and the 
second in 6 ways, all of which ways are equally likely to hap- 
pen, then the probability that the first will happen is expressed 
mathematically by the fraction a/(a + 4), while the probability 
that the second will happen is 6/(@+6). Such events are said 
to be mutually exclusive. Denote their probabilities by p and 
g respectively. Then there result 


a b 
p-37p (age 
the last equation following from the first two and being the 


mathematical expression for the certainty that one of the two. 
events must happen. 


Thus, to illustrate, in tossing a coin it must give “ head ” or 
“tail”;@=6=1,andg=g=1/2. Again, if an urn contain 
a@=5 white and 6 = 8 black balls, the probability of drawing 


* Jevons, Principles of Science, New York, 1874, p. 217. 
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a white ball in one trial is = 5/13 and that of drawing a 
black one g = 8/13. 

Similarly, if there are several mutually exclusive events 
which can occur in a, b,c... ways respectively, their probabil- 
ities 2, 9,7... are given by 


a ae c 
en. Fie...’ adele 
(2) 
P2t+t¢et7+...=1. 


For example, if an urn contain a= 4 white, d= 5 black, 
and ¢ = 6 red balls, the probabilities of drawing a white, black, 
and red ball at a single trial are p= 4/15, g = 5/15, and 
r = 6/15, respectively. 

Formulas (1) and (2) may be applied to a wide variety of 
cases, but it must suffice here to give only a few such. As a 
first illustration, consider the probability of drawing at random 
a number of three figures from the entire list of numbers which 
can be formed from the first seven digits. A glance at the 
table of Art. 1 shows that the symbols of formula (1) have in 
this case the values a= 210, and a+6= 13699. Hence 
6 = 13489, and p = 210/13699; that is, the probability in ques- 
tion is about 1/65. 

Secondly, what is the probability of holding in a hand of 
whist all the cards of one suit? Formula (1) of Art. 3 shows 
‘that the number of different hands of thirteen cards each which 
may be formed from a pack of fifty-two cards is 


Bg. 61<S0%...40 
egigecuiy oe 

and the probability required is the reciprocal of this number. 
The probability against this event is, therefore, very nearly 
unity. 

Thirdly, consider the probabilities presented by the case of 
an urn containing 4 white, 5 black, and 6 red balls, from which 
at a single trial three balls are to be drawn. Evidently the 
triad of balls drawn may be all white, all black, all red, partly 


white and black, partly white and red, partly black and red, or 
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one each of the white, black, and red. There are thus seven 
different probabilities to be taken into account. The theory 
of combinations shows (see equation (1), Art. 3) that the total 
number of 


White triads na; seaoaee ==) Ave 
Id 

Black triads ae 2 Sa IO re 

Red triads ae 2 at = 20=c 
: ie & 8.7 - 

White and black triads = = =( 441) =o 
: : 10. 

White and red triads ees Dae 4+20)== 96 =e 

li. 10.9 
Black and red triads = 5 = (10-F 20)=ags Fo) 
White, black, and red triads = 4.5.6 = 130, 


Sum = 455 
The total number of these triads is 455, and is, as it should 
be, the number of combinations in groups of three each of the 
whole number of balls. Hence formulas (2) give the seven 
different probabilities which follow, using the initial letters 
_w, 6, r to indicate the colors represented in a triad: 


ope ec aie P= 4/455, 
CO ThE) 6 bbb I= 10/455, 
py Y=: “20/455, 
“6 6% aezyh or wbhb 5 70/455; 
“ 6 6 any OF WIT t= 96/455, 
« «bby or brr u = 135/455, 
“6&6 yhy v = 120/455. 


Prob. 7. When three dice are thrown together, what is the prob- 
ability that the throw will be greater than g? 

Prob. 8. Write down a literal formula for the saghabiiaes of the 
several possible triads considered in the above question of the balls, 
supposing the numbers of white, black, and red balls to be Z, m,n, 
respectively. 
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ART. 5. PROBABILITY OF CONCURRENT EVENTS. 


If the probabilities of two independent events are f, and 
Py» respectively, the probability of their concurrence in any 
single instance is ~,f, Thus, suppose there are two urns 
U, and U,, the first of which contains a, white and 4, black 
balls, and the second a, white and 3d, black balls. Then the 
probability of drawing a white ball from UV, is ~,=a,/(a, + 4,), 
while that of drawing a white ball from UV, is p, = a,/(a, + 6,). 
The total number of different pairs of balls which can be formed 
from the entire number of balls is (a, ++ 4,)(a, + 4,). Of these 
pairs a,a, are favorable to the concurrence of white in simul- 
taneous or successive drawings from the two urns. Hence the 
probability of a concurrence of 

white with white = a,a,/(a, + 6,)\(a,+ 4,), 

white with black = (a,4,-+ a,0,)/(a, + 4,)(a, + 8,), 

black with black = 0,6,/(a, + 6,)(a, + 4,), 
and the sum of these is unity, as required by equations (2) of 
Article 4. 

In general, if 2,, 2,, ~, . . . denote the probabilities of several 
independent events, and P denote the probability of theit 
concurrence, 

Pee DED DS. (1} 

To illustrate this formula, suppose there is required the 
probability of getting three aces with three dice thrown simul- 
taneously. In this case ~, =p, = ~, = 1/6 and 
Fea (1/6)° = 1/216. 

Similarly, if two dice are thrown simultaneously the proba- 
bility that the sum of the numbers shown will be 11 is 2/36; 
and the probability that this sum 11 will appear in two succes- 
sive throws of the same pair of dice is 4/36.36. 

The probability that the alternatives of a series of events 
will concur is evidently given bv 

O=49.9.9,---=(1-AXI—P/)MI—-p,)--- (2) 

Thus, in the case of the three dice mentioned above, the 

probability that each will show something other than an ace is 
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2, = 9% = 9%, = 5/6, and the probability that they will concur in 
this is Q = 125/216. 

Many cases of interest occur for the application of (1) and 
(2). One of the most important of these is furnished by suc 
cessive trials of the same event. Consider, for example, what 
may happen in z trials of an event for which the probability 
is f and against which the probability is g. The probability 
that the event will occur every time is p”. The probability that 
the event will occur (#7 — 1) times in succession and then fail is 
p'>'qg. But if the order of occurrence is disregarded this last 
combination may arrive in different ways; so that the prob. 
ability that the event will occur (z— 1) times and fail once is 
up*~*'g. Similarly, the probability that the event will happer 
{x — 2) times and fail twice is $u(z — 1) p"~*g’; etc. That is, 
the probabilities. of the several possible occurrences are given by 
the corresponding terms in the development of (f + gq)”. 


By the same reasoning used to get equations (2) of Art. 
3 it may be shown that the maximum term in the expansion 
of (p+ )” is that in which the exponent m, say, of g is 
the whole number lying between (z+ 1)g — 1 and (z+ 1)g. 
In other words, the most probable result in » trials is the 
occurrence of the event (~—m) times and its failure m 
times. When ~ is large this means that the most probable of 
all possible results is that in which the event occurs  — ng 
= n(I — g) = up times and fails mgtimes. Thus, if the event 


be that of throwing an ace with a single die the most probable: 


of the possible results in 600 throws is that of 100 aces and 
500 failures. 

Since g” is the probability that the event will fail every time 
in 2 trials, the probability that it will occur at least once in x 
trials is 1 —g*. Calling this probability 7,* 


r=1-— Gi =a (3) 
If vy in this equation be replaced by 1/2, the corresponding 
value of is the number of trials essential to render the 


* See Poisson’s Probabilité des Jugements, pp. 40, 41. 


ee  ———— 


. PROBABILITY OF CONCURRENT EVENTS, 21 


chances even that the event whose probability is p will occur 
at least once. Thus, in this case, the value of » is given by 

os log 2 

TA ees 
This shows, for example, if the event be the throwing of double 
sixes with two dice, for which p = 1/36, that the chances are 
even (ry = 1/2) that in 25 throws (z = 24.614 by the formula) 
double sixes will appear at least once. 


Equation (3) shows that however small may be, so long as 
it is finite, 2 may be taken so large as to make 7 approach in- 
definitely near to unity; that is, 2 may be so large as to render 
it practically certain that the event will occur at least once. 


When ~ is large 


mwn—t) , u(z—i)(n—2), 
a ia ae 


{1 —p)"* = 1r— np + ——— ae Wied eag 


= e~"* approximately. 


sree 


Thus an approximate value of 7 is 


Tie Fang A, log ¢ = 0.4342495. (4) 
This formula gives, for example, for the probability of drawing 
the ace of spades from a pack of fifty-two cards at least once in 
104 trials y = 1 — e~* = 0.865, while the exact formula (3) 
gives 0.867. 

Similarly, the probability of the occurrence of the event at 
least ¢ times in # trials will be given by the sum of the terms 
of (p+ ¢)” from p” up to that in #’g"~‘ inclusive. This proba- 
bility must be carefully distinguished from the probability that 


the event will occur ¢ times onlyin the x trials, the latter being 
expressed by the single term in f’g"~‘. 


Prob. 9. Compare the probability of holding exactly four aces in 
five hands of whist with the probability cf holding at least four aces 
in the same number of hands. 

Prob. 10. What is the probability of an event if the chances are 
even that it occurs at least once in a million trials ? See equation (4). 
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ArT. 6. BERNOULLI’S THEOREM. 


Denote the exponents of f and g in the maximum term of 

(p + 9)" by wand m respectively, and denote this term by 7. 
Then 
TH n(n — 1)(# — 2)... (M+ DD pgm a n} 


um! ~ pm! 


pg. (1) 


As shown in Art. 5, # in this formula is the greatest whole 
number in (z-+1)f, and m the greatest whole number in 
(x -+1)g; so that when x is large, wand m are sensibly equal 
to mp and vq respectively. 


The direct calculation of 7 by (1) is impracticable when z 
is large. To overcome this difficulty the following expression 
is used : * 


n\aante* Yam + oot i+... (2) 


12” 
log ¢e = 0.4342495, log 27 =0.7981799. 


This expression approaches n"e-" ¥27u as a limit with the 
increase of #, and in this approximate form is known as Stir- 
ling’s theorem. Although a rude approximation to mz! for 
small values of # this theorem suffices in nearly all cases 
wherein such probabilities as ZT are desired. Making use of 
the theorem in (1) it becomes 


(inp ne 


V2anpgq (3) 


That this formula affords a fair approximation even when 
m is small is seen from the case of a die thrown 12 times. The 
probability that any particular face will appear in one throw is 
p = 1/6, whence g = 5/6; and the most probable result in 12 
throws is that in which the particular face appears twice and 
fails to appear ten times. The probability of this result com- 
puted from (3) is 0.309, while the exact formula (1) gives 0.296. 


The probability that the event will occur a number of times 


* This expression is due to Laplace, Théorie Analytique des Probabilités, 
See also De Morgan’s Calculus, pp. 600-604, 
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comprised between (4 — a) and (4+ a) in x trials is evidently 
expressed by the sum of the terms in (f+ g)" for which the 
exponent of g has the specified range of values. Calling this 
probability A, putting 

h=nmp+u, and m=ng—4, 


and using Stirling’s theorem (which implies that ” is a large 
number),* 


I u \~ “b+ u \~9-%) 
raz ginl+ a)" 2) 
Y2mnpg Lh up : nq y 


very nearly; and the summation is with respect to u from 
“u—=—atou=-+a. But expansion shows that the natural 
logarithm of the product of the two binomial factors in this 
equation is approximately — u’*/2npg. Hence 


I 


V2nnpq 


—42 
R= oO pMe « 


and, since z is supposed large, this may be replaced by a definite 
integral, putting 


dz=1/V2npq, and 2° = u'/2npg. 


Thus 
+a/ Vinpy a/ Y2npq 
K= real aa £ Bm ef cee (4) 
—a/ 9p : 


This equation expresses the theorem of James Bernoulli, 
given in his Ars Conjectandi, published in 1713. 

The value of the right-hand member of (4) varies, as it 
should, between o and 1, and approaches the latter limit rap- 
idly as z increases. Thus, writing for brevity 


2 
i=-— = de. 
"ead eed 
0 


* See Bertrand, Calcul des Probabilités, Paris, 1889, for an extended discus- 
sion of the questions considered in this Article. 
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the following table shows the march of the integral : 


z I F 4 ig Zz Hs 
0.00 0.000 0.75 0.711 1.50 0.966 

+25 .276 I.00 .843 1.75 .987 

-50 .520 1.25 -923 2.00 -995 


To illustrate the use of (4), suppose there is required the 
probability that in 6000 throws of a die the ace will appear a 
number of times which shall be greater than 1/6 K 6000 — I0 
and less than 1/6 X 6000-++ IO, or a number of times lying 
between ggo and 1010. In this case a = 10, 2 = 6000, p = 1/6, 
g= 5/6. Thus, a/V2upg = 10/¥2.6000.1/6.5/6 = 0.245. 
Hence, by (4) and the table, R = 0.27. 


Prob. 11, If the ratio of males to females at birth is 105 to 100, 
what is the probability that in the next 10,000 births the number of 
males will fall within two per cent of the most probable number ? 


Prob. 12. If the chance is even for head and tail in tossing a 
coin, what is the probability that in a million throws the difference 
between heads and tails will exceed 1500? 


ART. 7. INVERSE PROBABILITIES.* 


If an observed event can be attributed to any one of several 
causes, what is the probability that any particular one of these 
causes produced the event? To put the question in a concrete 
form, suppose a white ball has been drawn from one of two 
urns, U, containing 3 white and 5 black balls, and U, contain- 
ing 2 white and 4 black balls; and that the probability in favor 
of each urn is required. If U, is as likely to have been 
chosen as U,, the probability that U, was chosen is 1/2. After 
such choice the probability of drawing a white ball from U, is 
3/8. Before drawing, therefore, the probability of getting a 
white ball from U, was 1/2 X 3/8 = 3/16, by Art. 5. Similarly, 
before drawing the probability of getting a white ball from U, 
was 1/2 X 2/6=1/6. These probabilities will remain un- 
changed if the number of balls in either urn be increased or 


* See Poisson, Probabilité des Jugements, pp. 81-83. 
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diminished so long as the ratio of white to black balls is kept 
constant. Make these numbers the same for the two urns. 
Thus let the first contain g white and 15 black, and the second 
8 white and 16 black; whence the above probabilities may be 
written 1/2 X 9/24 and 1/2 X 8/24. It is now seen that there 
are (9+ 8) cases favorable to the production of a white ball, 
each of which has the same antecedent probability, namely, 1/2. 
Since the fact that a white ball was drawn excludes considera- 
tion of the black balls, the probability that the white ball came 
from U, is 9/17 and that it came from U, is 8/17; and the sum 
of these is unity, as it should be. 

To generalize this result, let there be # causes, C,, C,,... C,,. 
Denote their direct probabilities by 7,,9,,...9,,; their antecedent 
probabilities by 7,, 7,,...7,; and their resultant probabilities 
on the supposition of separate existence by /,,2,,...- Dm: 

That is, 
PRET ee Rey og epost? (1) 

Let D be the common denominator of the right-hand mem- 
bers in (1), and denote the corresponding numerators of the 
several fractions by s,, 5,,...5,. Then 

Pee get Pita SSD), ya Py, Se Le 
and it is seen that there are in all (s,-+5,-+...5,,) equally 
possible cases, and that of these s, are favorable to C,, s, to 
Peeuaence, it F., P,... £.denote the probabilities, of 
the several causes on the supposition of their coexistence, 
I =s/(,+5,+- + Sm) =P,/(P, + Ps ++ + + Pm) 

Thus in general 

P.=p,/2p, P,= b./ 2p, ++» Pa =Pn/ 2p. (2) 

To illustrate the meaning of these formulas by the above 

concrete case of the urns it suffices to observe that 
fone. | geet 4/5 and 7, == 1/2, 
for U,, 97,= 1/3 and r,= 1/2; 
whence ~,= 3/16, 2,= 1/6, ~,+2.= 17/48; 
and 2 ewOt, Fo = 8/47. 


As a second illustration, suppose it is known that a white 
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ball has been drawn from an urn which originally contained 
balls, some of them being black, if all are not white. What is 
the probability that the urn contained exactly white balls? 
The facts are consistent with m different and equally probable 
hypotheses (or causes), namely, that there were 1 white and 
(m—1) black balls, 2 white and (m#— 2) black balls, etc. 


Hence‘in (1), g, = 4,'=... = 1, and 
Poi lm, Py 2], cs fn SUL eo hag Se 
Thus =p = (1/2)(m-+ 1), 
2n 
and Pes pf mim + 1)" 


This shows, as it evidently should, that 2 = m is the most 
probable number of white balls in the urn. The probability 
for this number is P,, = 2/(m -++ 1), which reduces, as it ought, 
to1form=1. — 

Formulas (1) and (2) may also be applied to the problem of 
estimating the probability of the occurrence of an event from 
the concurrent testimony of several witnesses, X,, X,,... 
Denote the probabilities that the witnesses tell the truth by 
4, x, ... Then, supposing them to testify independently, 
the probability that they will concur in the truth concerning 
the event is 7,7,...; while the probability that they will con- 
cur in the only other alternative, falsehood, is (I—+,)(I—~,) ... 
The two alternatives are equally possible. Hence by equations 


(1) and (2) 


Aya Rityien =O EM — aoe 

P= Pe A 

+) Ft ees 

‘pl (I— +4,)(1 —%,) ... (3) 


Bg a = ee 
P, being the probability for and P, that against the event. 

To illustrate (3), if the chances are 3 to 1 that_X, tells the 
truth and 5 to 1 that_X, tells the truth, +, = 3/4, x, = 5/6, and 
P,= 15/16; or, the chances are 15 to 1 that an event occurred 
if they agree in asserting that it did.* 


* For some interesting applications of equations (~) see note E of Appendix 
to the Ninth Bridgewater Treatise by Charles Babbage (London, 1838). 
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It is of theoretical interest to observe that if x,, 7,,...in 
(3) are each greater than 1/2, P, approaches unity as the 
number of witnesses is indefinitely increaséd. 


Prob, 13. The groups of numbers of one figure each, two figures 
each, three figures each, etc., which it is possible to form from the 
nine digits 1, 2,...g9 are printed on cards and placed severally in 
nine similar urns. What is the probability that the number 777 will 
be drawn in a single trial by a person unaware of the contents of 
the urns? 


Prob. 14. How many witnesses whose credibilities are each 3/4 
are essential to make P, = 0.999 in equation (3) ? 


ArT. 8. PROBABILITIES OF FUTURE EVENTS. 


Equations (2) of Art. 7 may be written in the following 
manner: 


RCE Namah mie ©) 
If 2,, 2.) + + «Pm are found by observation, P,, P,,...P,, will ex- 
press the probabilities of the corresponding causes or their 
effects. When, as in the case of most physical facts, the num- 
ber of causes and events is indefinitely great, the value of any 
p or P in (1) becomes indefinitely small, and the value of >p 
must be expressed by means of a definite integral. Let x de- 
note the probability of any particular cause, or of the event to 
which it gives rise. Then, supposing this and all the other 
causes mutually exclusive, (I — x) will be the probability 
against the event. Now suppose it has been observed that in 
(m +n) cases the event in question has occurred m times and 
failed x times. The probability of such a concurrence is, by 
Art. 5,cx"(1 — x)", where c is a constant. Since + is unknown, 
it may be assumed to have any value within the limits 0 and 1; 
and all such values are a priori equally possible. Put 

x= cx™"(1 — 2)". 

Then evidently the probability that + will fall within any as- 
signed possible limits a and 4 is expressed by the fraction 


J ss / of ea: 
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so that the probability of any particular x is given by 


a1 — 4)"dx 


yy a"(1 — x)"ax 


This may be regarded as the antecedent probability of the 
cause or event in question. 


ot (2) 


What then is the probability that in the next (7 + s) trials 
the event will occur 7 times and fail s times, if no regard is had 
of the order of occurrence? If # were known, the answer 
would be by Arts, 2 and 5 

(r+s)1_ : 

teers a’"(1 — +) (3) 
But since x is restricted only by the condition (2), the required 
probability will be found by taking the product of (2) and (3) 
and integrating throughout the range of x. | Thus, calling the 
required probability Q, 


amr — x)"Vde 


Gest e sia (4) 
vi a1 — x)'dx 


The definite integrals which appear here are known as Gamma 
functions. They are discussed in all of the higher treatises on 
the Integral Calculus. Applying the rules derived in such 
treatises there results * 


oat Inte 49) (m4n41)! 
~ rls!mini(mtutrts4ti)! ° (5) 
If regard is had to the order of occurrence of the event; 


that is, if the probability required is that of the event happen- 
ing 7 times in succession and then failing s times in succession, 


* It is a remarkable fact that formula (5) is true without restriction as to 
values of m, 2,7, s. The formula may be established by elementary considera- 
tions, as was done by Prevost and Lhuilier, 1795. See Todhunter’s History of 
he Theory of Probability, pp. 453-457. 
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the factor (r+ s)!/r!s! in (3), (4), (5) must be replaced by 
unity. 

To illustrate these formulas, suppose first that the event 
has happened m times and failed no times. What is the prob» 
ability that it will occur at the next trial? In this case (4) 
gives 


Q = furde] f xax = (m+ 1)/(m-+ 2). 


When m is large this probability is nearly unity. Thus, the 
sun has risen without failure a great number of times m; the 
probability that it will rise to-morrow is 


I 7 ee I 2 
(1+ =)(r +2) Se cae Goentin We Wert 
which is practically 1. 


Secondly, suppose an urn contains white and black balls in 
an unknown ratio. If in ten trials 7 white and 3 black balls 
are drawn, what is the probability that in the next five trials 
2 white and 3 black balls will be drawn? The application of 
(5) supposes the ratio of the white and black balls in the urn 
to remain constant. This will follow if the balls are replaced 
after each drawing, or if the number of balls in the urn is sup- 
posed infinite. The data give 


m7, n= 3, r=z2, sS= 3, 
mtr=9 ut+s=6, rt+s=5, m+n+I1=-I1, 
mtnutrt+s+i1= 16. 
Thus by (5) 


s!o!6!r1! 
O=21317131161~ 39/9 


Suppose there are two mutually exclusive events, the first 
of which has happened m times and the second # times in 
m-—+-n trials. What is the probability that the chance of the 
occurrence of the first exceeds 1/2? The answer to this ques- 
tion is given directly by equation (2) by integrating the nume- 
rator between the specified limits of x. That is, 
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if: a™(1 — x)"ax 
Poe ©) 


i (1 — ae 


Thus, if # = 1 and x =0, P= 3/4; or the odds are three to 
one that the event is more likely to happen than not. Simi- 
larly, if the event has occurred m times in succession, 


Ved ——90 Week fr) oa A 


which approaches unity rapidly with increase of x. 


ART. 9. THEORY OF ERRORS. 


The theory of errors may be defined as that branch of math- 
ematics which is concerned, first, with the expression of the re- 
sultant effect. of one or more sources of error to which com- 
puted and observed quantities are subject ; and, secondly, with 
the determination of the relation between the magnitude of 
an error and the probability of its occurrence. In the case of 
computed quantities which depend on numerical data, such as 
tables of logarithms, trigonometric functions, etc., it is usually 
possible to ascertain the actual values of the resultant errors. 
In the case of observed quantities, on the other hand, it is not 
generally possible to evaluate the resultant actual error, since 
the actual errors of observation are usually unknown. In either 
case, however, it is always possible to write down a symbolical 
expression which will show how different sources of error enter 
and affect the aggregate error; and the statement of such an 
expression is of fundamental importance in the theory of errors. 


To fix the ideas, suppose a quantity Q to be a function of 


several independent quantities x, y, z...; that is, 

= HEI, face 
and let it be required to determine the error in Q due to errors 
in z, y, 2... Denote such errors by 40, Ar, Ay, Az... 


Then, supposing the errors so small that their squares, prod- 
ucts, and higher powers may be neglected, Taylor’s series gives 
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49 = 2 ae 4 Cay 4 a4. (1) 
_ This equation may be said to express the resultant actual error 
of the function in terms of the component actual errors, since 
the actual value of JQ is known when the actual errors of 
x,y,2...are known. It should be carefully noted that the 
quantities 4, y, ... are supposed subject to errors which are 
independent of one another. The discovery of the independent 
sources of error is sometimes a matter of difficulty, and in general 
requires close attention on the part of the student if he would 
avoid blunders and misconceptions. Every investigator in work 
of precision should have a clear notion of the error-equation of 
the type (1) appertaining to his work; for it is thus only that 
he can distinguish between the important and unimportant 
sources of error. 


Prob. 15. Write out the error-equation in accordance with (r) 
for the function Q = xyz + x* log (y/2). 

Prob. 16. In a plane triangle a/6 = sin A/sin B. Find the error 
in @ due to errors in 4, A, and &. 

Prob. 17. Suppose in place of the data of problem 16 that the 
angles used in computation are given by the following equations : 
A=4A,+4(180°—4,— 8,—C,), B= B, + 4(180° — 4,—B, - C)), 
where 4,, &,, C, are observed values. What then is Ja? 

Prob. 18. If w denote the weight of a body and ¢ the radius of 
the earth, show that for small changes in altitude, 47 /w=—2J4r/r; 
whence, if a precision of one part in 500 000 000 is attainable in com- 
paring two nearly equal masses, the effect of a difference in altitude 
of one centimeter in the scale-pans of a balance will be noticeable.* 


ART. 10. LAWS OF ERROR. 


A law of error is a function which expresses the relative 
frequency of occurrence of errors in terms of their magnitudes. 
Thus, using the customary notation, let e€ denote the magni- 


* This problem arose with the International Bureau of Weights and Measures, 
whose work of intercomparison of the Prototype Kilogrammes attained a pre- 
cision indicated by a probable error of 1/500 000 oooth part of a kilogramme. 
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tude of any error in a system of possible errors. Then the law 
of such system may be expressed by an equation of the form 
I = (6): (1) 
Representing € as abscissa and y as ordinate, this equation 
gives a curve called the curve of frequency, the nature of which, 
as is evident, depends on the form of the function ¢ This | 
equation gives the relative frequency of occurrence of errors in 
the system; so that if € is continuous the probability of the 
occurrence of any particular error is expressed by yde = o(e)de; 
which is infinitesimal, as it plainly should be, since in any con- 
tinuous system the number of different values of e is infinite. 


Consider the simplest form of ¢(e), namely, that in which 
~(€) =c,aconstant. This form of ¢(e€) obtains in the case of 
the errors of ‘tabular logarithms, natural trigonometric func- 
tions, etc. In this case all errors between minus a half-unit 
and plus a half-unit of the last tabular place are equally likely 
to occur. Suppose, to cover the class of cases to which that 
just cited belongs, all errors between the limits —aand +a 
are equally likely to occur. The probability of any individual 
error will then be ¢(€)de = cde, and the sum of all such prob- 
abilities, by equation (2), Art. 4, must be unity. That is, 


saa 4a 
J ede =e fide = 1. (2) 


This gives ¢ = 1/2a, or by (1) y=1/2a. The curve of fre- 

Q quency in this case is shown in the figure, 

C D AB being the axis of € and OQ that of y. 

It is evident from this diagram that if the 

errors of the system be considered with 

respect to magnitude only, half of them 

should be greater and half less than @/2. 

= o gs This is easily found to be so in the case of 
tabular logarithms, etc. 

As a second illustration of (1), suppose y and e connected 


by the relation y = ¢ Va* — e’, where a is the radius of acircle, 
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¢ a constant, and € may have any value between — a and + a. 
Then the condition 


cl deVa'—eé=1 


gives c= 2/(a’z). In this, as in the preceding case, (+ e) = 
g(— ¢), the meaning of which is that positive and negative 
errors of the same magnitude are equally likely'to occur. It 
will be noticed, however, that in the latter case small errors 
have a much higher probability than those near the limit a, 
while in the former case all errors have the same probability. 


In general, when e€ is continuous P(e) must satisfy the condi- 
tion J ode = I, the limits being such as to cover the entire 


range of values of e. The cases most commonly met with are 
those in which ¢(€) is an even function, or those in which 
f(+ e€) = ¢(—«). In such cases, if + a denote the limiting 
value of €, 


fede = 2 f ede —a5 (3) 


ART. 11. TYPICAL ERRORS OF A SYSTEM. 


\ 


Certain typical errors of a system have received special 
designations and are of constant use in the literature of the 
theory of errors. These special errors are the probable error, 
the mean error, and the average error. The first is that error 
of the system of errors which is as likely to be exceeded as 
not ; the second is the square root of the mean of the squares 
of all the errors; and the third is the mean of all the errors 
regardless of their signs. Confining attention to systems in 
which positive and negative errors of the same magnitude 
are equally probable, these typical errors are defined mathe- 
matically as follows. Let 

€, = the probable error, 
€,,== the mean error, 
€, = the average error. 
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Then, observing (2) of Art. 10, 
—& +ep 


aE P(e\de = -f sede ae P(ej)de = =i P(e\de = 


a / ‘Hleede, €, = 2 Ve p(ejede. 


The student should seek to avoid the very common misap- 
prehension of the meaning of the probable error. It is not 
“the most probable error,” nor “the most probable value of 
the actual error”; but it is that error which, disregarding signs, 
~vould occupy the middle place if all the errors of the system 
were arranged in order of magnitude. A few illustrations will 
suffice to fix the ideas as to the typical errors. Thus, take the 
simple case wherein P(e) =c= 1/2a, which applies to tabular 
jogarithms, etc. Equations (1) give at once 


(1) 


&—=—+ =a, eo ape i + <a. 

For the case of tabular values, a= 0.5 in units of the last 

tabular place. Hence for such values 
6é= 10:25, €& = 20.20, €,=" 2025. 

Prob. 19. Find the typical errors for the cases in which the law 
of error is P(e) = Va? — &’, P(e) =e tafe), P(e)=c cos*(ze/2a); 
¢ being a constant to be determined in each case and € having any 
value between — a and + a. 


ArT. 12. LAWS OF RESULTANT ERROR. 


When several independent sources of error conspire to pro- 
duce a resultant error, as specified by equation (1) of Art. 9, 
there is presented the problem of determining the law of the 
resultant error by means of the laws of the component errors. 
The algebraic statement of this problem is obtained as follows 
for the case of continuous errors: 

In the equation (1), Art. 9, write for brevity 


30Q 30 
e= JQ, e, = —4z, =—=4 gevie 5 
Q = € 5 Ly . 
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and let the laws of error of e, €,, €,,... be denoted by ¢(e),. 
p,(€,), %,(€,)... Then the value of e¢ is given by 
e=e,te,+... (1) 
The probabilities of the occurrence of any particular values. 
of €, €,,...are given by ¢,(e¢,)de,, ¢,(e,)de,,...; and the 
probability of their concurrence is the probability of the cor- 
responding value of e. But since this value may arise in an 
infinite number of ways through the variations of Gee Ges net 
over their ranges, the probability of e, or (e)de, will be 
expressed by the integral of ¢,(e,)de,,(e,)de, ... subject to 
the restriction (1). This latter gives e, = e— e, —e,..., and 
de, = de for the multiple integration with respect to e,, e,,..- 
Hence there results 


Ple)de = de [| o,(€ —€, —&,—. . .),(e,)de,..., 
or 


o(6) = [dle—2— 2, —.. Jola)de, fdleddes--- (2) 


It is readily seen that this formula will increase rapidly im 
complexity with the number of independent sources of error.* 
For some of the most important practical applications, how- 
ever, it suffices to limit equation (2) to the case of two inde- 
pendent sources of error, each of constant probability within 
assigned limits. Thus, to consider this case, let e, vary over 
the range —a to +a, and e, vary over the range — 6 to + 6. 
Then by equation (2), Art. 10, 

¢,(€,) = 1/(2a),  o,(€,) = 1/(20). 


Hence equation (2) becomes 


I 2 
P(e) = real deé,. 


In evaluating this integral e€, must not surpass + 6 and 
€, = € — €, must not surpass + a. Assuming a > 4, the limits 
of the integral for any value of e=e,-+ ¢, lying between 
— (a+ 6) and — (a — 4) are —dand+(e+a). This fact is 

* The reader desirous of pursuing this phase of the subject should consult 


Bessel’s Untersuchungen ueber die Wahrscheinlichkeit der Beobachtungsfehler; 
Abhandlungen von Bessel (Leipzig, 1876), Vol. II. 
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made plain by a numerical example. For instance, suppose 
a@=sandé=3. Then —(a+6)=—8and —(a— 4)= — 
Take e = — 6,a number intermediate to — 8 and — 3. Then 
the following are the possible integer values of e, and e, which 
will produce e = — 6: 


€ €, Ee limits of e, 
—6=—5—1, —1=+(e+4) 
acai. east 


ae tS ets Nn ys aatar - 
Similarly, the limits of e, for values of e lying between 
— (a — 6) and + (a— d) are — 6 and + 34; and the limits of 
e, for values of € between + (a— 0) and + (a+ 6) are + (€—a) 
and +46. Hence 
e+a 


HO i5 | as for —(a+6) < e< — (a—6), 


-b 


+b 


5 
HO= 75 | d= yan for =(a— 8) ex ee @) 


+b 
H9=75 J de =— SESH for + (a0) <e< (049). 


a 


Thus it appears that in this case the graph of the resultant 
law of error is represented by the upper base and the two sides 
of a trapezoid, the lower base being the 
axis of e and the line joining the middle 
points of the bases being the axis of $(e). 
(See the first figure in Art. 13.) The prop- 
erties of (3), including the determination 
of the limits, are also illustrated by the 
adjacent trapezoid of numerals arranged 
to represent the case wherein a = 0.5 and 
6=0.3. The vertical scale, or that for P(e), does not, how- 
ever, conform exactly to that for e. 


IIOI! 
IITIOIII 
ELPIORTLS 
Bert Orr Bik 
ELITIIOLLIIITE 
LPELITEIORI PELE 
PIGTELSIOLEEELSE 
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Prob. 20. Prove that the values of #(€) as given by equation (3) 
satisfy the condition specified in equation (3), Art. 10. 

Prob, 21. Examine equations (3) for the cases wherein a = Jand 
4= 0; and interpret for the latter case the first and last of (3). 

Prob. 22. Find from (3), and (1) of Art. rz, the probable error of 
the sum of two tabular logarithms. 


ART. 138. ERRORS OF INTERPOLATED VALUES. 


Case I.—One of the most instructive cases to which formulas 
{3) of Art. 12 are applicable is that of interpolated logarithms, 
trigonometric functions, etc., dependent on first differences. 
Thus, suppose that v, and v, are two tabular logarithms, and 
that it is required to get a value v lying ¢ tenths of the interval 
from v, towards v, Evidently 

v=v,+(v, —v,) t= (1 — Zu, + 4, ; 
and hence if ¢, ¢,, e, denote the actual errors of v, v,, v, , re- 
spectively, 
e=(1 — Ze,+ ¢e,. (1) 

It is to be carefully noted here that ¢ as given by (1) re- 
quires the retention in v of at least one decimal place be- 
yondthe last tabular place. For example, let v = log (24373) 
from a 5-place table. Then v,= 4.38686, v, = 4.38703, 
v, — Vv, = +0.00017, ¢ = 0.3, and v = 4.38691.1. Likewise, as 
found from a 7-place table, e, = — 0.45, e, = + 0.37in units of 
the fifth place; and hence by (1)¢= — 0.20. That is, the 
actual error of v = 4.38691.1 is = 0.20, and this is verified by 
reference to a 7-place table. 

The reader is also cautioned against mistaking the species 
of interpolated values here considered for the species common- 
ly used by computers, namely, that in which the interpolated 
value is rounded to the nearest unit of the last tabular place. 
The latter species is discussed under Case II below. 


Confining attention now to the class of errors specified by 
equation (1), there result in the notation of the preceding 


article 
é,=(1—Z)e, €& =e, and e=e—e,+6,; 


and since ¢, and ¢, each vary continuously between the limits 
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+o.5 of a unit of the last tabular place, a and @ in equations 
(3) of that article have the values 
a=0.5(f —%), 6=0.54 
Hence the law of error of the interpolated values is ex- 
pressed as follows: 


g(€) = a ae for values of e betw. —o.5 and —(0.5—2), | 


ah ane” for values of € betw. —(0.5 —Z) and +-(0.5 —4), | (2) 


5 O5—€ 


= GM for values of € betw. +(0.5—2) and -+-0.5. 


The graph of ¢(e) for ¢= 1/3 is shown by the trapezoid 
AB, BC, CD in the figure on page 40, Evidently the equa- 
tions (2) are in general represented by a trapezoid, which degen- 
erates to an isosceles triangle when ¢ = 1/2. 

The probable, mean, and average errors of an interpolated 
value of the kind in question are readily found from (2), and 
from equations (1) of Art. 11, to be 


ge, = (1/4\(1 — 2) fot’ .0 <<a 
= 1/2 — (1/2)W2at(1 — 2) for 1/3 <¢ < 2/3, 
met OE forez {3° <7 <1, 
Lat = 277 | 
soe" oe or — OF ae r (3) 
1 — (1 — 2#)° 


= a rage fOF Ore aye: 


1 — (2¢ — 1)* 
Yagi — ae 
It is thus seen that the probable error of the interpolated 
value here considered decreases from 0.25 to 0.15 of a unit of 
the last tabular place as ¢ increases from 0 to 0.5. Hence such ° 
values are more precise than tabular values; and the computer 
who desires to secure the highest attainable precision with a 
given table of logarithms should retain one additional figure 
beyond the last tabular place in interpolated values. 


for-2/ 2 <oete | 
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Case II.—Recurring to the equation v = v, + ¢(v,— v,) for an 
interpolated value v in terms of two consecutive tabular values 
v, and v,, it will be observed that if the quantity ¢(v,— v,) is 
rounded to the nearest unit of the last tabular place, a new error 
is introduced. For example, if v,=log 1633 = 3.21299, and 
v, = log 1634 = 3.21325 from a 5-place table, v, — v,= +26 
units of the last tabular place; and if ¢= 1/3, ¢(v, — v,) = 83; 
so that by the method of interpolation in question there results 
V = 3.21299 + 9 = 3.21308. Now the actual errors of v, and 
v, are, as found froma 7-place table, — 0.38 and +-0.21 in units 
of the fifth place. Hence the actual error of uv is by equation 
(1), ?X —0.38-+4 x +0.21 —4= — 0.52, as is shown di- 
rectly by a 7-place table. 

It appears, then, that in this case the error-equation cor- 
responding to (1) is 

e=(1—Ze,+ te, +8,, (4) 
wherein e, and e, are the same as in (1) and ¢, is the actual error 
that comes from rounding ¢(v, — v,) to the nearest unit of the 
last tabular place. 

The error ¢,, however, differs radically in kind from e, and 
e, The two latter are continuous, that is, they may each have 
any value, between the limits — 0.5 and +-0.5; while ¢, is dis- 
continuous, being limited to a finite number of values depend- 
ent on the interpolating factor 7. Thus, for ¢= 1/2 the only 
possible values of e, areo-+ 1/2, and — 1/2; likewise for ¢ = 
1/3, the only possible values of e, are 0, + 1/3, and — 1/3. It 
is also clear that the maximum value of e, which is constant and 
equal to 1/2 for (1), is variable for (4) in a manner dependent 
on?#. For example, in (4), 

The maximum of e=1/2-+1/2=1, fort=1/2, 
“ He =1/2 41/3 = 5/6, “ t= 1/3, 
«“ “ “emaf2--31/2—1, “ i= 1/4; 
+ s “eat /e4-2/5-= 0/10" > t= 1/5. 

The determination of the law of error for this case presents 
some novelty, since it is essential to combine the continuous 
errors (I — fe, and ¢e, with the discontinuous error ¢,. The 
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simplest mode of attacking the problem seems to be the fol- 
lowing quasi-geometrical one. In the notation of Arts. 12 and 
13, put in (4) e= 6 (I— Ze, = €,, te, = €,, and e, =e, Then 

e=(e,+.6¢,) +6, (5) 

The law of error for (€, + €,) is given by equation (2) for any 
value of 4. Hence for a given value of ¢ there will be as many 
expressions of @(e) as there are different values of €,. The 
graphs of $(e) will all be of the same form but will be differently 
placed with reference to the axis of P(e). Thus, if ¢= 1/3 the 
values of e€, are — 1/3, O, and 
+ 1/3, and these are equally 
likely to occur. For e,=0 the 
graph is given directly by (2), 
and is the trapezoid ABCD 
\ symmetrical with respect to OQ. 
\ For ¢,=-—1/3 the graph is 
\ abQd, of the same form as 
‘+ ABCD but shifted to the left 
by the amount of e,=— 1/3. 
Similarly, the graph for the case 
of ¢, = + 1/3 is a’Qd’d’, and is produced by shifting ABCD to 
the right by an amount equal to + 1/3. 

Now, since the three systems of errors for this case are 
equally likely to occur, they may be combined into one system 
by simple addition of the corresponding element areas of the 
several graphs. Inspection of the diagram shows* that the 
resultant law of error is expressed by 

P(e) = (1/4)(5 + Oe) for — 5/6 <e<— 1/6, 
=I for—1/6<e<+1f6, (6) 
= (1/4(5—6e) for + 1/6 <e< 45/6 


This is represented by a trapezoid whose lower base is 10/6, 
upper base 2/6, and altitude I. 


/ 
/ \ 
a [Ds & 


* Sum the three areas and divide by 3 to make resultant area = I, as 
required by equation (3), Art. 10. 
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As a second illustration, consider equation (5) for the case 
¢= 1/2. In this case e, must be either o or 1/2, the sign of 
which latter is arbitrary. For e, = 0, equations (2) give 

P(e) =2+4€ for —1/2<€<0, 
=2-—4e for ENE (7) 
This function is represented by the isosceles triangle AQE 
whose altitude OQ is twice the base AZ. 

Similarly #(e) for e, = + 1/2 would 


be represented by the triangle AQE dis- 
placed to the right a distance 1/2; and 
if the two systems for €,=0 and e, = 
+ 1/2 be combined into one system, 
their resultant law of error is evidently 
P(e)=1+2e for —1/2<€e<0, 

=I for o< e<+ 1/2, (8) 

= 2—2e for+1/2<e<1; 4 7 5 5 
the graph of which is ABCD. On the 
other hand, if the errors inthis combined system be considered 
with respect to magnitude only, the law of error is 

g(e) = 2(1—e) for o<e<tI, (9) 

the graph of which is OQD. 

The student should observe that (6), (7), (8), and (9) satisfy 


the condition J welde = 1 if the integration embraces the 


whole range of e. 

The determination of the general form of ¢(e) in terms of 
the interpolating factor ¢ for the present case presents some 
difficulties, and there does not appear to be any published solu- 
tion of this problem.* The results arising from one phase of 
the problem have been given, however, by the author in the 
Annals of Mathematics,t and may be here stated without 
proof. The phase in question is that wherein ¢ is of the form 
1/n, n being any positive integer less than twice the greatest 


* The author explained a general method of solution in a paper read at the 
summer meeting of the American Mathematical Society, August, 1895. 


+ Vol. II, pp. 54-59. 
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tabular difference of the table to which the formulas are ap- 
plied. For this restricted form of ¢ the possible maximum 
value of € as given by equation (5) is, in units of the last 
tabular place, (2% — 1)/n for ” odd and 1 for ” even. 


The possible values of €, of equation (5) are 


I 2 t= 
re) _ SS Tac for z odd 
is + + + 2n ‘ 
oe ae cee wits f 
fe) ~ ae - for 2 even. 
: n° n’ ay 2 2 


An important fact with regard to the error 1/2 for z even 
is that its sign is arbitrary, or is not fixed by the computation 
as is the case with all the other errors. However, the com- 
puter’s rule, which makes the rounded last figure of an inter- 
polated value even when half a unit is to be disposed of, will, 
in the long-run, make this error as often plus as minus. 


The laws of error which result are then as follows: 


For x odd. 
g(€) = 1 for e between — 1/2” and + 1/2n, 
n (2n—1 
d(e) = or ( + e) for ebetw. F 1/2” and F(2u—1)/2n., 
For 7 even. 
—2 
P(€) = eI + e) for e between o and F 1/n, 
n Geose - 
oa eae + €} for € betw. F 1/m and + (zm — 1)/n, 
= oe poy + €) for e between F (#—1/m) and $1. 


By means of these formulas and (1) of Art. 11 the probable, 
mean, and average errors for any value of x can be readily 
found. The following table contains the results of sucha com- 
putation for values of z ranging from 1 to 10. The maximum 
actual error for each value of z is also added. The verifica- 
tion of the tabular quantities will afford a useful exercise to the 
student. 


"sae = 
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TypIcAL Errors OF INTERPOLATED LOGARITHMS, ETC. 


Interpotating | Probable Error. Mean Error. Average Error. 


tor. Maxi 
pa ah “s : Em co Kesani yess 
I 0.250 0.289 0.250 1/2 
1/2 292 408 333 I 
1/3 -256 +347 +287 5/6 
1/4 .276 ae tea| i313 j I 
1/5 .268 +370 -303 9/10 
1/6 -277 +385 315 I 
9G 274 .380 PSU 13/14 
1/8 * 3270 -389 -318 I 
1/9 -278 .386 -316 17/18 
1/10 .281 +392 +320 I 


When the interpolating factor ¢ has the more general form 
m/n, wherein m and m are integers with no common factor, the 
possible values of €, are the same as fort=1/z. But equa- 
tions (3) of Art. 12 are not the same for ¢= m/n as fort = 1/n, 
and hence for the more general form of ¢, @(€) assumes a new 
type which is somewhat more complex than that discussed 
above. ‘The limits of this work render it impossible to extend 
the investigation to these more complex forms of #(€). It may 
suffice, therefore, to give a single instance of such a function, 
namely, that for which ¢= 2/5. For this case 

P(€) = 1 for e between 0 and F 1/10, 
(5/6)(13/10 + e) for e between F 1/10 and F 3/10, 

= (5/3)(4/5 + €) for e between F 3/10 and F 7/I0, 
= (5/6)(9/10+ €) for e between F 7/10 and + 9/Io. 

The graph of the right-hand half of 4) 8 
this function is shown in the accompany- (el 
ing diagram, the whole graph being 
symmetrical with respect to OA, or the 
axis of (€). 

Attention may be called to the strik- 
ing resemblance of this graph to that of 
the law of error of least squares. o 


E 


Prob. 23. Show from equations (3) that €,, varies from 1/12 
= 0.29 —, for = 0, to 1/¥24 =0.20+, for ¢= 0.5; and that €, 
varies from 0.25 to 1/6 for the same limits. 
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Prob. 24. Show that the probable, mean, and average errors 
for the case of ¢= 2/5 cited above (p. 43) are + 0.261, + 0.251, 
and + 0.290, respectively. 


ArT. 14. STATISTICAL TEST OF THEORY. 

A statistical test of the theory developed in Art. 13 may 
be readily drawn from any considerable number of actual er- 
rors of interpolated values dependent on the same interpolating — 
factor. The application of such a test, if carried out fully by 
the student, will go far also towards fixing clear notions as to 
the meaning of the critical errors. 

Consider first the case in which an interpolated value falls 
midway between two consecutive values, and suppose this 
interpolated value retains two additional figures beyond the 
last tabular place. Then by equations (2), Art. 13, the law of 
error of this interpolated value is 

P(e) = 2-++ 4e for € between — 0.5 and o 
= 2 — 4e for e between 0 and-+-0.5. 


Hence by equation (1) of Art. 11, or equation (3) of Art. 12, the 
probable error in this system of errors is 4 — (4)V¥2=0.15. 
It follows, therefore, that in any large number of actual errors 
of this system, half should be less and half greater than 0.15. 
Similarly, of.the whole number of such errors the percentage 


falling between the values 0.0 and 0.2 should be 

+o.2 +o.2 

J 6d = of (2 — 4€)de = 0.64; 
that is, sixty-four per cent of the errors in question should be 
less numerically than 0.2. 

To afford a more detailed comparison in this case, the act- 
ual errors of five hundred interpolated values from a 5-place 
table have been computed by means of a 7-place table. The 
arguments used were the following numbers: 20005, 20035, 
20065, 20105, 20135, etc., in the same order to 36635. The 
actual and theoretical percentages of the whole number of 
errors falling between the limits 0.0 and 0.1, 0.1 and 0.2, etc., 
are shown in the tabular form following: 
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Panis wi Rcrore: Actual Theoretical 

Percentage. Percentage. 
SARA OL is ceigtetes sane s 6 ssc cess ese 3362 36 
RMN wikls aisles <'e)s ie a isincn s > ae x “302 28 
REMINDED ctw a nial Salpti'e 40 0 -0,'k 50 Se «- I9.0 20 
Me Oe eens as cet aoc aide a+ 7132 12 
DA MAITG Os), 2 xis stern actin oh onite aiats a ae ree og Be 4 
CEVA OAS 15 is rs ose hues one pra 51.4 50 


The agreement shown here between the actual and theoretical 
percentages is quite close, the maximum discrepancy being 2.8 
and the average I.5 per cent. 

Secondly, consider the case of interpolated mid-values of the 
species treated under Case II of Art. 13. The law of error for 
this case is given by the single equation (9) of Art. 13, namely, 
P(e) = 2(1 — e), no regard being paid to the signs of the errors. 
The probable error is then found from 

&p 
af — e)de = 3, 
whence «, = 1 —4$¥2=0.29. Similarly, the percentage of 
the whole number of errors which may be expected to lie, for 
example, between 0.0 and 0.2 in this system is 


0.9 
2 f(s — e)de = 0.36. 
° 
Using the same five hundred interpolated values cited 
above, but rounding them to the nearest unit of the last tabu- 
lar place and computing their actual errors by means of a 7-place 
table, the following comparison is afforded : 


AA Actual Theoretical 
Limits of Errors. Percentage. Percentage. 
0.0 and 0.2...... Meee SORTER es Loe * 36 
ee BN O.4 cave ese Pas aie u.awin's 2 27.8 28 
0.4 and 0.6....cecsecsece PA cee tls 20 
G.O and 0.8.6 5..00: Peuguake tale vak aad 12 
SNIUISE Sey ages bien gcaves ces 5.6 4 


0.0 and 0.29....++. Sidemiates names 49.8 50 
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The agreement shown here between the actual and theoretical 
percentages is somewhat closer than in the preceding case, the 
maximum discrepancy being only 1.6 and the average only 0.6 
per cent. 

Finally, the following data derived from one thousand act- 
ual errors may be cited. The errors of one hundred inter- 
polated values rounded to the nearest unit of the last tabular 
place were computed * for each of the interpolating factors 
0.1, 0.2,...0.9. The averages of these several groups of act- 
ual errors are given along with the corresponding theoretical 
errors in the parallel columns below: 


Interpolating Actual Theoretical 
Factor. _ Average Error. Average Error, 
eR OP Pe Pee Three hy ee Gets 0.320 
Se dew wk gang Oe Lene ees ee ee oe 0.303 
OA a adiaiieds isso de eee Pe ee aes 0.304 
OAva Lnyen be otisn Pie vEs Centar ae 0.290 
oe PRERCER Tee A lee 0.333 
CED slate igen ni oictajh wa(e pedir turns as Ola} 0.290 
WF es ds mah oon os Waok wale ie One® 0.304 
OSvisi vee es en? ees eet eS) 0.303 
Ye i oauces er le aaceyes é valetnya? eee, 0.320 


The average discrepancy between the actual and theoret- 
ical values shown here is 0.017. It is, perhaps, somewhat 
smaller than should be expected, since the computation of the 
actual errors to three places of decimals is hardly warranted 
by the assumption of dependence on first differences only. 


The average of the whole number of actual errors in this 
case is 0.308, which agrees to the same number of decimals 
with the average of the theoretical errors. ¢ 


* By Prof. H. A. Howe. See Annals of Mathematics, Vol. III, p. 74. 
The theoretical averages were furnished to Prof. Howe by the author. 

+ The reader who is acquainted with the elements of the method of least 
squares will find it instructive to apply that method to equation (1), Art. 13, 
and derive the probable error of ¢. This is frequently done without reserve by 


STATISTICAL TEST OF THEORY. AT 


Prob. 25. Apply formulas (3) of Art. 12 to the case of the sum: 
or difference of two tabular logarithms and derive the correspond- 
ing values of the probable, mean, and average errors. The graph 
of #(€) is in this case an isosceles triangle whose base, or axis of €, 
is 2, and whose altitude, or axis of #(€), is 1. 


those familiar with least squares. Thus, the probable error of ¢ or ¢2 being 
0.25, the probable error of ¢ is found to be 


0.25 Vi — 2t+ 27, 


This varies between 0.25 for¢=o0 and 0.18 for¢=4; while the true value of 
the probable error, as shown by equations (3), Art. 13, varies from 0.25 to 0.15 
for the same values of ¢. It is, indeed, remarkable that the method of least 
squares, which admits infinite values for the actual errors ¢; and é, should give 
so close an approximate formula as the above for the probable error of e. 

Similarly, one accustomed to the method of least squares would be inclined 
to apply it to equation (4), Art. 13, to determine the probable error of e. The 
natural blunder in this case is to consider ¢, ¢2, and ¢s independent, and ¢s like 
¢, and é continuous betweer the limits 0.0 and 0.5; and to assign a probable 
error of 0.25 to each. In tis manner the value 


0.25 V2(r — ¢+ 2%) 


is derived. But this is absurd, since it gives 0.25 V2 instead of 0.25 for ¢= 0. 
The formula fails then to give even approximate results except for values of ¢ 
near 0.5. 


ee 
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